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ABSTRACT 

We present the magnetic and electric field structures as well as the currents and charge 
densities of pulsar magnetospheres which do not obey the ideal condition, E • B = 0. Since the 
acceleration of particles and the production of radiation requires the presence of an electric field 
component parallel to the magnetic field, E||, the structure of non-Ideal pulsar magnetospheres 
is intimately related to the production of pulsar radiation. Therefore, knowledge of the structure 
of non-Ideal pulsar magnetospheres is important because their comparison (including models 
for the production of radiation) with observations will delineate the physics and the parameters 
underlying the pulsar radiation problem. We implement a variety of prescriptions that support 
nonzero values for Ey and explore their effects on the structure of the resulting magnetospheres. 
We produce families of solutions that span the entire range between the vacuum and the (ideal) 
Force-Free Electrodynamic solutions. We also compute the amount of dissipation as a fraction 
of the Poynting flux for pulsars of different angles between the rotation and magnetic axes and 
conclude that this is at most 20-40% (depending on the non-ideal prescription) in the aligned 
rotator and 10% in the perpendicular one. We present also the limiting solutions with the 
property J = pc and discuss their possible implication on the determination of the "on/off" 
states of the intermittent pulsars. Finally, we find that solutions with values of J greater than 
those needed to null Ey locally produce oscillations, potentially observable in the data. 

Subject headings: Magnetohydrodynamics (MHD)-Methods: numerical-pulsars: general-Gamma rays: 
stars 



1. Introduction 

Pulsars are extraordinary objects powered by 
magnetic fields of order 10^^ G anchored onto neu- 
tron stars rotating with periods ^ 10^'^ — 10 s. 
These fields mediate the conversion of their rota- 
tional energy into MHD winds and the accelera- 
tion of particles to energies sufficiently high to pro- 
duce GeV photons. Their electromagnetic emis- 
sion is quite complex and ranges from the radio to 
the multi-GeV 7— ray regime. 

Pulsar radiation has been the subject of many 
studies since their discovery; however, despite 
such a multi-decade effort by a large number of 
researchers its details are not fully understood. 



There are several reasons for this. On the ob- 
servational side, the sensitivity of the 7— ray in- 
struments (the band where most of the pulsar ra- 
diation power is emitted) has become adequate 
for the determination of a significant number of 
pulsar spectra only recently. Since its launch in 
2008, the Fermi Gamma-Ray Space Telescope has 
revolutionized pulsar physics by detecting nearly 
90 7— ray pulsars so far, with a variety of well- 
measured light curves. On the theoretical side, 
the structure of the pulsar magnetosphere, even 
for the simplest axisymmetric case, remained un- 
certain for a long time. The modeling of pulsar 
magnetospheres has seen rapid advancements only 
very recently. 
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The original treatment of the no n-axisymmetric 
rotating stellar dipole in vacuum (jPeutschl 119551 ) 
established the presence of an electric field com- 
ponent, E|j3, parallel t o the magnetic field on the 
star's surface. Then Goldreich fc Julian 
suggested that this component of the electric field 
would pull charges out of the pulsar, opening up 
magnetic field lines that cross the Light Cylinder 
(hereafter LC), to produce an MHD wind with cur- 
rents flowing out of and into the neutron star polar 
cap. The charge density needed to achieve this is 
known as the Goldreich - Julian density pc J ; an as- 
sociated current density Jqj — cpGi is obtained if 
we allow these charges to move outward at nearly 
the speed of light. 

Following the work of lGoldreich fc Julian (Il969l) 
a number of attempts were made to produce the 
structure of the pulsar magne tosphere . Thus, 



Scharlemann fc Wagoned (|1973I ) were able to re- 



duce the structure of the axisymmetric magneto- 
sphere to a single equation for the poloidal mag- 
netic flux, the so-called pulsar equation. However 
an exact solution of this equation that would pro- 
duce the magnetic field structure across the LC 
was missing for almost 3 decades (for a review 
of the problems and the various attempts to pro- 
duce a consistent pulsar magnetosphere , even in 
the axisymmetric case, see Michell 1982 ). In fact, 
the inability of models to determine the magne- 
tospheric structure across the LC led people to 
speculate that its structure is not smooth across 
this surface and assumed that the observed radia- 
tion is the res ult of magnetic field disc ontinuities 
across the LC (|Mestel fc Shibatalll99l ). 

The structure of the axisymmetric pulsar 
magnetosphere within Force Free Electrodynam- 
ics (hereafter FFE) (i.e. assuming that the 
only surviving component of the electric field 
is that perpendicular to the magnetic field, 
or equivalently that E ■ B = ) was given by 
IContopoulos. Kazanas. fc Fendtl ( 1999t ) (hereafter 
CKF) who used an iterative procedure to deter- 
mine the current distribution on the LC so that 
the magnetic field lines should cross this surface 
smoothly. Besides allowing a smooth transition of 
the magnetic field lines across the LC, this solution 
also resolved the issue of the current closure in pul- 



sar magnetospheres: The current starts near the 
polar axis and closes, for R > i?Lc (here Rhc is the 
radius of the LC), mainly along an equatorial cur- 
rent sheet (allowing also a small amount of return 
current at a finite height above it). For R < JJlc 
the flow is mainly along the separatrix of the open 
and closed field lines. The axisymmetric solution 
has since been confirmed and further studied b'v 



nas smce oeen connrmea ana lu rtner stua i ea by 
several others (fC ruzinov 2O O5I: iTimokhinI hOOd 



Komissarovl bOOalMcKinnevI I2OO6I) . in particular 



with respect to its properties near the Y-point 
(the int ersection of the last closed field line and 
the LC: IUzdensk"^l2003l ). 

The first non-axisymmetric (3D), oblique rota- 
tor m agnetosphere was presented by ,Spitkov skvl 
(|2006l ). who used a time-dependent numerical code 
to advance the magnetic and electric fields under 
FFE conditions to steady state. These simulations 
confirmed the general picture of current closure 
established by the CKF solution and produced a 
structure very similar to that of CKF in the ax- 
isymmetric case^_Simi]A£_simu^^ 
formed bv iKalapo tharak os &: ContopoulosI ( 20091 ) 
who, using a Perfectly Matched Layer (hereafter 
PML) at the outer boundary of their computa- 
tional domain, were able to follow these simula- 
tions over many stellar rotation periods. In gen- 
eral, the 3D magnetosphere, just like the axisym- 
metric one, consists of regions of closed and open 
field lines with a large scale electric current circuit 
established along open magnetic field lines. In the 
3D case, the current sheet needed for the global 
current closure is in fact u ndulating, a,s fore seen 
in the kinematic solution of Bogovalov ( 19991 ) . As 
shown by recent simulations, the undulating cur- 
rent sheet structure is stable at least to distances 
of IORlc and for several stellar rotation period s 



^ E|| = (E ■ B)B/B^, where E, B are the electric and mag 
netic field vectors. 



(jKalapotharakos. Contopoulos. fc Kazanasll201ir) . 

The FFE solutions presented so far, by con- 
struction, do not allow any electric field compo- 
nent parallel to the magnetic field B (E • B = 0); 
as such they do not provide the possibility of par- 
ticle acceleration and the production of radia- 
tion, in disagreement with observation. Nonethe- 
less, models of pulsar radiation have been con- 
structed assuming either a vacuum dipole ge- 
ometry, or the field line geometry of the FFE 
magnetosphere. Polar cap particle accelera- 
tion, 7— ray emission and the formation of pair 
cascades were discussed over many years as 
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a means of acco unting for the observed high- 
energy radiation ('Ruder man fc Sutherland! 1975t 
Daughertv fc Hard ing 199^. However, 



measure- 



ments by Fermi o f the cutoff shape of the Vela 



pulsar spectrum (|Abdo et alj l2009bl ) has ruled 



out the super-exponential shape of the spectrum 
predicted by polar cap models due to the atten- 
uation of pair production. Models placing the 
origin of the high-energy emission in the outer 
magnetosphere or beyond are now universally 
favored. Slot gap model s ( Muslimov fc Harding! 
2004 : Harding et al.!!2008h . with acceleration and 



emission in a narrow layer along the last open 
field line from the neutron star sur face to near the 
LC, and outer gap niqdels (jll omani fc Yadigaroglu 



1995!: !Romani!ll99l !Takata et al.l!2007[) produce 



light curves and spectr a in first-order agreement 
with the Fermi data (Romaiii fc Watters 2010t 
!Venter. Harding, fc Guillemot! !2009! ) . All of the 
above models assume a retarded vacuum magnetic 
field. There are also kinematic models of pulsar 
emission derived by injecting photons along the 
magnet ic field lines of either the FFE magneto- 
sphere ( Bai fc Spitkovskv 2010b|) or those of th e 
vacuum field geometry ( Bai fc Spitkovskv 2010al ). 
or along the electric current near the equatorial 
current sheet where the ratio J/ pc — » 1 in FFE 



(jContopoulos fc KalapotharakosI l2010l ). These 
models also produce light curves in broad agree- 
ment with observations. However, fits of Fermi 
data light curves with slot gap and outer gap mod- 
els in FFE geometry are less favor able than those 



for t he vacuum dipole geometry (jHarding et al 



20111 ) 



While FFE magnetospheres do not allow for 
particle acceleration, they do, however, provide 
a global magnetospheric structure consistent with 
the boundary conditions on the neutron star sur- 
face (a perfect conductor), a smooth transition 
through the LC and the establishment of a large 
scale MHD wind. They also determine both the 
currents and the sign of the charges that flow in 
the magnetosphere. However, they do not pro- 
vide any information about the particle produc- 
tion that might be necessary to support the under- 
lying charge conservation and electric current con- 
tinuity. This was apparent already i n the C KF so- 
lution and was also noted in .Michel! (|l982r) : There 
are regions, most notably the intersections of the 
null charge surface with open magnetic field lines 



and the region near the neutral (Y) point of the LC 
where the FFE charge density changes sign along 
an open field line (along which particles are sup- 
posed to flow freely) . The numerical FFE solution 
has no problem in these regions because it can sim- 
ply supply at will the charges necessary for current 
closure and for consistency with the global bound- 
ary conditions. This, however, requires sources of 
electron-positron pairs at the proper regions and 
of the proper strength. Even though the FFE 
treatment of the pulsar magnetosphere involves no 
such processes, these must be taken into account 
self-consistently. Viewed this way, pulsar radia- 
tion is the outcome of the "tension" between the 
boundary conditions imposed by the global mag- 
netospheric structure and the particle production 
necessary for charge conservation and electric cur- 
rent continuity. 

In Figs. IT^,b,c we present the structure of the 
FFE magnetosphere on the poloidal plane fi — fl 
for inchnation angles a = (0°, 45°, 90°). These fig- 
ures depict the poloidal magnetic field lines and in 
color the value of J/pcH The vertical dashed lines 
denote the position of the LC. In the axisymmet- 
ric case the value of J/pc is less than one in most 
space and exceeds this value only in the null sur- 
faces and the separatrix between open and closed 
field lines. The value 1 of this ratio is significant 
because it denotes the maximum current density 
for carriers of single sign charge. Values of this 
ratio greater than 1, by necessity denote the pres- 
ence of both positive and negative charges moving 
in opposite directions. For a ~ 0, J/pc < 1 in most 
of the polar cap region, however, as the value of 
inclination a increases so does the size of regions 
with J/pc > 1; in the perpendicular rotator we 
see clearly that J / pc> 1 over the entire polar cap 
region. 

With the above considerations in mind, the 
question that immediately arises is how would 
the structure of the pulsar magnetosphere and 
the consequent production of radiation be modi- 
fied, if the ideal condition, namely E • B = 0, were 
dropped. This is the subject of the present work. 
A successful magnetospheric configuration (solu- 
tion) should be able to address the following ques- 
tions: a) Where does particle acceleration take 



^Hereafter, wherever we use the ratio J/pc we consider its 
absolute value. 
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place? b) How close are these processes to a 
steady state? c) Is the ensuing particle acceler- 
ation and radiation emission consistent with the 
Fermi pulsar observations? 

Unfortunately, while the ideal condition is 
unique, its modifications are not. The evolution 
equations require a relation between the current 
density J and the fields, i.e. a form of Ohra's law , 
a relation which is currently lacking. iMeieii (j2004r ) 



proposed a generalization of Ohm's law in plasmas 
taking into account many effects (e.g. pressure, 
inertia. Hall). However, in the present case we 
are in need of a macroscopic expression for Ohm's 
law which would be also everywhere compatible 
with the microphysics self-consistently, a proce- 
dure that might involve Particle in Cell (PIC) 
techniques, well beyond the scope of this paper. 
For this reason, the present work has a largely ex- 
ploratory character. A number of simple prescrip- 
tions are invoked which lead to magnetospheric 
structures with values of E|| ^ 0. The global elec- 
tric current and electric field structures are then 
examined with emphasis on the modifications of 
the corresponding FFE configuration. We are 
presently interested only in the broader aspects 
of this problem, deferring more detailed calcula- 
tions to future work. In particular, the general 
compatibility of the E|| magnitude with the val- 
ues needed to produce (through pair production) 
self-consistent charge and current densities from 
the plasma microphysics will not be discussed in 
the present paper. The contents of the present 
paper are as follows: In §2 we discuss the differ- 
ent modifications in the expression for the electric 
current that we apply in order to obtain non-Ideal 
solutions. In §3 we present our results and the 
properties of the new solutions. In §4 we give 
particular emphasis to specialized solutions that 
obey certain global conditions. Finally, in §5 we 
present our conclusions. 

2. Non-Ideal Prescriptions 



Maxwell equations 

— = -cV X E 
dt 

aE 

— = cV X B - AttJ 
dt 

under ideal force-free conditions 



(1) 
(2) 



E-B = 0, 



pE + - J X B 
c 



0, 



where p = V • E/(47r). The evolution of these 
equations in time requires in addition an expres- 
sion for the current density J as a function of E 
and B. This is given by 



E X B 



B-VxB-E-VxE 



' = '^^"-4^ ^ 

()Gruzinovlll9"99l) . The second term in Eq. ^ en- 
sures that the condition E ± B (the ideal condi- 
tion) is preserved during the time evolution. How- 
ever, this term introduces numerical instabilities 
since i t involves spatial d erivatives of both fields. 
Thus, ISpitkovskvl (|2006l ) implemented a scheme 



that evolves the fields considering only the first 
term of Eq. ([3]) , and at the end of each time step 
"kills" any developed electric field component par- 
allel to the magnetic field (Em). In the case of 



non-Ideal status, the E 



condition does not 



apply anymore and electric field components par- 
allel to the magnetic field are allowed. However, 
while in FFE the prescription that determines the 
value of E|| is unique (it is equal to zero), there is 
no unique prescription in the non-Ideal case. Be- 
low we enumerate the non-Ideal prescriptions that 
we consider in the present study: 
(A) The above implementation of the ideal con- 
dition, hints at an easy generalization that leads 
to non-Ideal solutions: One can evolve Eqs. ([T]) & 
([2]), using only the first term of the FFE current 
density (Eq. [3]) , and at each time step keep only a 
certain fraction b of the E|| developed during this 
time instead of forcing it to be zero. In general 
the portion b of the remaining E|| can be either 
the same everywhere or variable (locally) depend- 

p, J). As b goes 



■ nom to 1 the corresponding solution goes from 



In t he FFE descript ion of pulsal magneto- . , , / 

, 1 ^ . ■ : — I ] „"" „ k ,1 .. , ; , — — — iiig on some other quantity (e.g 

spheres lopitkovskvl ( 20061 1 and lKalapotharakos fc Contopo ulor 
~i ' ' ' ^—^^^^^^^^^—^^^^^^^—^ j^^.yj 

( 20091 ) solved numerically the time dependent ^ t ii • • r 

' ' I* i< E to vacuum. In this case an expression for the 

electric current density is not given a priori, and 

J can be obtained indirectly from the expression 



J = -LfeVxB-^ 

47r V dt 



(4) 
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(B) Another way of controlling E|| is to introduce 
a finite conductivity a. In this case we replace the 
second term in Eq. ([3]) by crE|| and the current 
density reads 



cp- 



E X B 



52 



crE 



II 



(5) 



Note that Eq. ([5]) is related to but is not quite 
equivalent to Ohm's law whic h is defined in 
the frame of the fluid. Others (Lvutikovl 2003 : 



Li. Spitkovskv. fc Tchekhovskovl2oTlF implemented 
a different version closely related to Ohm's law. 
The problem is that the frame of the fluid is not 
well defined in our problem. Moreover, it is not a 
priori known neither that a single form of Ohm's 
law is applicable in the entire magnetosphere, nor 
that the conductivity corresponding to each such 
form of Ohm's law should be constant (e.g. in 
regions where pair production takes place). This 
is why we chose a simpler expression, which turns 
out to yield results qualitatively very similar to 
other formulations based on Ohm's law. As is the 
case in prescription (A) with the parameter 6, here 
too, as a ranges from cr = to cr ^ oo we obtain 
a spectrum of solutions from the vacuum to the 
FFE, respectively. We note that even though the 
vacuum and FFE solutions are unambiguously 
defined limiting cases, there is an infinite num- 
ber of paths connecting them according to how a 
depends on the local physical parameters of the 
problem {p,E,B etc.). 

(C) It ha s been argued (|Lvubarskvll996t l iGruzinovl 
2007t lLvubarskv..2008i) that regions corresponding 
to space- like currents {J/pc > 1) are dissipative 
due to instabilities related to counter streaming 
charge flowsH Moreover, the space-like current 
regions should trace t he pair produc tion areas 
(i.e. dissipative areas). ICruzinov ( 2007 ) proposed 
a covariant formulation for the current density 
that introduces dissipation only in the space-like 
regions while the time-like regions remain non- 
dissipative. The current density expression in the 
so called Strong Field Electrodynamics (hereafter 
SEE) reads 

^ _ CpE X B + (c2p2 + ^2^2^2)1/2(5^3 + SqE) 



(6) 



^Note that space-like currents (J/pc > 1) are formed exclu- 
sively by counter streaming charge flows. 



where 



- El = B^ - E^ BoSo = E . B, £;o > (7) 
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El 



Bl 



El 



(8) 



and CT is a function of Eq , Bq with dimensions of 
conductivity (inverse time). The expression for 
J (Eq. ini) produces either a space-like (u > 0) 
or a null [a = 0) current. Note that this treat- 
ment precludes time-like currents. Such currents 
do exist but they are only effectively time-like: 
the current J, along with E\\ , fluctuates contin- 
uously, so that while J is locally greater than (or 
equal to) pc its average value is less than that, 
i.e. (J) < pc. This behavior is captured by our 
code and is manifest by the continuously changing 
direction of the parallel component of the electric 
fleld along and against the magnetic held direction 
()Gruzinovil20"08ll201ll) . 

In prescriptions (A) and (B) above, care must 
be taken numerically so that the resulting value 
of E± (defined as |E - E|||) be less than B. Vi- 
olation of this condition usually happens in the 
current sheet where the value of B goes to zero. 
We note that the presence of the Eq term in the 
denominator of the SEE presc ription (Eg. El) and 
of the prescription of J used bv lLi et al.l (|201lh in- 
deed prevents the denominator from going to zero 
and alleviates this issue in the current sheet re- 
gions. However, alleviation of this problem does 
not guarantee that all treatments that include this 
Eq term converge to the same solution or that 
even produce the same dissipation rates in the cur- 
rent sheet. Thus, the solutions and in p articular 



the di ssipation rates of SEE and that of iLi et al 
(j201l[ ) can be quite different. 

All the simulations that are presented be- 
low have run in a cubic computational box 
20 times the size of the LC radius -Rlc, i-e 
[— 10i?Lc ■ ■ • 10i?Lc]'^- Outside this area we imple- 
mented a PML layer that allows us to follow the 
evolution of the magnetosphere for several stella r 
rotations (jKalapotharakos fc Contopoulos [iooi). 
The adopted grid cell size is 0.02i?LC, the stellar 
radius has been considered at — 0.3i?LC, and 
all simulations were run for 4 full stellar rotations 
to ensure that a steady state has been achieved. 
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3. Results 

We have run a series of simulations using pre- 
scription (B) that cover a wide range of a valuefl 
In Fig. [2^ we present the values of the Poynting 
flux L integrated over the surface of the star for the 
a = 0° (aligned; red), a = 45° (green) and a = 90° 
(perpendicular; blue) versus a in log-linear scale. 
The two solid horizontal line segments of the same 
color denote in each case (a = 0°, 45°, 90°) the cor- 
responding limiting values of the Poynting flux for 
the FFE and the vacuum solutions, respectively. 
The Poynting flux as a function of the inclina- 
tion angle a for the vacuum and the FFE solutions 
reads 



L 



Vacuum 



(9) 



(1 + si\? a) FFE 



where /z, f2 are the magnetic dipole moment and 
the rotational frequency of the star, respectively. 
This means that for a specific value of a the cor- 
responding value of the Poynting flux (spin-down 
rate) of each non- Ideal modification should lie be- 
tween the two values of Eq. Thus, Fig. [2^ 
shows that for the perpendicular rotator (and any 
case with a 7^ 0) the ratio of the Poynting flux over 
that of the vacuum reaches near 1 for low values 
of G (u/ri <i 0.01; for a = this is true only for 
cr — >■ 0). In Fig. [^b we plot for these simulations 
the corresponding dissipation power 



D 



J - E dF 



1 

47r 



ri<r<r2 



ri <r<r2 



(10) 

taking place within the volume bounded by the 
spheres ri = 0.3Rlc, i-e just above the star 
surface, and r2 = 2.5Rlc, versus a. We see 
that the energy losses due to dissipation exhibit 
a maximum that occurs at intermediate values of 
a {a/il ~ 1) for a = 90° and at higher values of 



''We have found from our simulations that the results of 
prescription (A) are very similar to those of prescription 
(B). We therefore restrict the discussion in the rest of this 
section to magnetospheres constructed using prescription 
(B). An explicit reference to results under prescription (A) 
is nonetheless made later on in this section in relation to 
results shown in Fig. 5. 



cr as the inclination angle a decreases. The dis- 
sipation power Ejj goes to as a goes either to- 
wards or 00 since both the vacuum (cr — 0) and 
FFE (cr — > 00) regimes are (by definition) dissipa- 
tionless. We note that despite the fact that the 
Poynting flux (on the stellar surface) varies signif- 
icantly with cr, the corresponding energy loss due 
to dissipation is limited and it never exceeds the 
value ~ 0.14/^^51^/0"^. The surface Poynting fiux 
reflects the spin down rate while the dissipative 
energy loss within the magnetosphere reflects the 
maximum power that can be released as radiation. 
Moreover, when Eu is measured as a fraction of 
the corresponding surface Poynting fiux it is al- 
ways less than 22% for a = 0°, 15% for a = 45° 
and 8% for a = 90°. These maximum fraction val- 
ues occur at lower a values than those that cor- 
respond to the maximum absolute values of Ejj 
(Fig.[lb,c). 

We note that all the values presented in Fig. [5] 
are the corresponding average values over 1 stellar 
period (4*^ stellar rotation). Moreover, we note 
that the values of Ed shown may be a little bit 
smaller than the true ones. This is because our 
simulations introduce some numerical dissipation 
even for the ideal-FFE solutions; for our resolu- 
tion, this produces a decrease in the Poynting flux, 
AL, which at r = 3 can reach ~ 10% of its value 
L on the stellar surface. This means that there is 
some energy loss that can not be measured by the 
J • E expression for the values provided by our so- 
lutions. Li et al. 2011 using higher resolution sim- 
ulations seem to obtain a corresponding decrease 
of « 5%. Due to this effect, the Poynting fiux 
reduction with r in the dissipative simulations is 
always larger than its true value. This means that 
the true Ed value is also a little higher than com- 
puted, since the energy reservoir is a little larger. 
Nevertheless, the error in Ed is less than AL/L. 
We ran also simulations with half grid cell size and 
simulations with the stellar radius set at 0A5Rlc- 
The Ed values in these cases were at most a few 
percent (mostly for small values of a) of the corre- 
sponding L values lower due to the finer descrip- 
tion of the stellar surface. 

In Fig. [3] we present the magnetospheric struc- 
tures that we obtain on the poloidal plane (/x, n) 
using prescription (B) for cr ~ 2417. Each row 
corresponds to the structure at a different incli- 
nation angle a = (0°,45°,90°) as indicated in 
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the figure. The first column shows the modu- 
lus of the poloidal current Jp in color scale to- 
gether with its streamlines. The second and third 
columns show, in color scale, the charge density 
p and E\\ respectively, together with the stream- 
lines of the poloidal magnetic field Bp. We ob- 
serve that for this high a value the global struc- 
ture s are similar to those of the FFE so lutions 
(see Contopoulos fc Kalapotharakos 20101 ). Even 
for the a = 0° case that has the most different 
value of the Poynting flux from that of the FFE 
solution, the current sheet both along the equator 
(outside the LC) and along the separatrices (inside 
the LC) does survive despite the effects of non- 
zero resistivity. However, we do observe magnetic 
field lines reconnecting gradually on the equato- 
rial current sheet beyond the LC. The main differ- 
ences between these solutions and the FFE ones 
are highlighted in the third column where we plot 
the parallel electric field E\\ developed during the 
evolution of the simulation. For a — Q° there is a 
significant parallel electric field component along 
the separatrices as well as over the polar capqj. 
The direction of the E|| follows the direction of the 
current J. Thus, E|| within the return current re- 
gion (the separatrices and a small part of the polar 
cap near them) points outward while in the rest of 
the polar cap region points inward. There is also 
a weaker parallel component E|| along the equato- 
rial current sheet (not clearly shown in the color 
scale employed). As a increases, the inward point- 
ing E|| component in the central part of the polar 
cap becomes gradually offset. During this topo- 
logical transformation the branch of the polar cap 
area corresponding to the return current becomes 
narrower while the other one becomes wider. At 
a = 90°, E|| is directed inwards in half of the polar 
cap and outwards in the other half. This topolog- 
ical behavior is similar to that of the poloidal cur- 
rent in the FFE solutions. For a ^ 0° there is a 
clearly visible component of E|| in the closed field 
lines area as well as along the equatorial current 
sheet outside the LC. 

Figure 0] shows the same plots as Fig. |3] and 
in the same color scale but for a ~ 1.517. Al- 
though the non-Ideal effects have been amplified, 
the global topological structure of the magneto- 



^We note that the concepts of the separatrix and the polar 
cap are less clearly in the case of non-Ideal solutions since 
magnetic field lines close even beyond the LC. 



sphere has not changed dramatically. We can still 
distinguish traces of the current sheet even in the 
a = 0° case, whose Poynting flux is substantially 
reduced from that of the FFE case. As noted 
above, the magnetic field lines close now beyond 
the LC, due to the finite resistivity that allows 
them to slip through the outflowing plasma. The 
parallel electric field component E|| is non-zero 
in the same regions as those seen in the higher 
(T— values; its topology also is not very different 
from that of these cases, but its maximum value 
increases with decreasing a. We note also that de- 
spite the fact that the dissipative energy losses cor- 
responding to the various values of a are quite dif- 
ferent (see Fig. [5^), the corresponding maximum 
values of are quite similar. 

We have also run many simulations implement- 
ing prescription (A) for various values of the frac- 
tion b of the non-zero parallel electric field com- 
ponent . We have covered this way the entire 
spectrum of solutions from the vacuum to the FFE 
one. The ensemble of solutions that connects these 
two limits seems to be similar to that of prescrip- 
tion (B). It seems that there are pairs of the values 
of the fraction h and the conductivity a that lead 
to very similar (qualitatively) solutions. In Fig. [5] 
we plot the poloidal magnetic field lines together 
with the parallel electric field component E\\ in 
color scale (similar to the third column of Figs. |3] 
and HI for h = 0.75 which has the same Poynting 
flux with the simulation of prescription (B) with 
a = 1.50. We see that the magnetospheric struc- 
tures of the two solutions are quite close (compare 
Fig. [5] with the third column of Fig. |4]). Also, 
ou r results are qualitatively very similar to those 
of Li et al. |2011) who, as noted, implemented a 
much more elaborate prescription for J based on 
Ohm's law. 

Figure |6] is similar to Fig. [2] but for prescrip- 
tion (C) - i.e. SEE, assuming the value of a to be 
constant. This prescription does not tend to the 
vacuum solution for (7 = but to a special solution 
with J/pc = 1. As in Figl2]all quantities presented 
in this figure correspond to average values over 1 
stellar period (4"^ rotation). The Poynting flux 
behaves accordingly, relaxing to L c± 0.55/i^ri^/c'^ 
for a = 0° (rather than zero), L ~ l.Qp?Q,'^/c? for 
a = 45° and L ~ 1.4^2^4/^3 ^j. a = 90° as cn> 0. 
These values are 55%, 65% and 70% of the cor- 
responding FFE values, respectively. For high a 
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values (fT oo) the corresponding solutions have 
a Poynting flux close to that of the FFE solutions. 
The dissipation power in this case is given by 

Ed= I J • E dF 



ri<r<r2 

and, as expected, is non-zero even for cr = 0. For 
(7 = we get the highest value for Ed, either 
in absolute sense or as a fraction of the Poynt- 
ing flux. These reach 45%, 25% and 12% for 
a — Q°,a = 45°, a — 90° respectively. These are 
higher than the highest values we got using pre- 
scription (B). Their increased magnitude is mostly 
due to the increase in the Ejj of the SFE solutions 
on the current sheet outside the LC. We observe 
also that the dissipation remains at high levels 
even for very high values of the conductivity and 
only the a — Q° case (which is the most dissipative 
one) shows a clear tendency of decreasing Ed with 
a. For (7 3> 1 the dissipation energy inside the LC 
decreases significantly while it remains almost con- 
stant outside the LC. We have tried higher resolu- 
tion simulations that allow higher values of a with- 
out any apparent decrease of the dissipation on the 
current sheet. This appears to confirm Gruzinov's 
claim that SFE allows for non-zero dissipation on 
the current sheet even for a ^ oo. 

Figure [7] depicts the same quantities as Figs.[3l 
|4]but for prescription (C) i.e. SFE and for a ~ 8f2. 
We observe that the magnetospheric structures in 
this case are very close to those of FFE (especially 
for a > 0°). The magnetic field lines open not 
much farther than the LC and the usual equato- 
rial current sheet develops beyond that. In the 
third column, which shows the parallel electric 
field, we see that in the effectively time-like re- 
gions near the star the color is "noisy" as a result 
of the constant change of its sign during the evo- 
lution. In the next section we are going to dis- 
cuss possible implications of these types of behav- 
ior in these regions. The regions with the "calm" 
(non noisy) electric fields correspond to space-like 
currents. We discern these electric fields mostly 
along the separatrix in the aligned rotator and 
along the undulating equatorial current sheet in 
the non-aligned rotator. We note that for a = 90° 
the entire polar cap region consists of such com- 



ponents. In Fig. |8] we present the results for the 
same prescription, i.e. (C) - SFE, but for a = 0. 
As we have mentioned in §2, these cases are sig- 
nificantly different from the vacuum solution. The 
magnetic field lines close well outside the LC while 
we are able to discern traces of the equatorial cur- 
rent sheet. We observe again the "noisy" regions 
near the star (this is the way that SFE deals with 
time-like currents, namely as fluctuating space-like 
ones). In the axisymmetric case (a = 0°) we see 
that the value of E'n as well as the area over which 
it occurs increase for the smaller values of a. For 
a > 0° the value that reaches near the polar 
cap is not significantly higher that that of Fig. [T] 
However, the new interesting feature in the oblique 
rotators is the development of high in specific 
regions of the closed field lines areas. These con- 
figurations do not imply high dissipation since the 
corresponding currents are very weak there. 

In order to reveal the on-average behavior of the 
effectively time-like regions we considered the av- 
erage values of the fields within one stellar period. 
In Fig. M we plot for the SFE (cr = 0) the cor- 
responding (average) parallel electric components 
together with the magnetic field lines (similar 
to the third columns in Figs. [7] and [5]). We see that 
the "noisy" behavior has disappeared. As we dis- 
cussed previously one can imagine the existence of 
hybrid solutions (e.g. non "noisy" behavior in the 
closed field line regions like that shown in Fig. |9]) 
and a "noisy" one in the open field line regions like 
that shown in Fig. [T) 

Finally, we note that for all the solutions pre- 
sented in this section, the regions corresponding 
to high values of J • E are traced mostly by the 
high values of J (e.g. near the current layers). 
However, as we observe in Figs. [3]|4l [7]|8] there are 
regions corresponding to high values of £'|| that do 
not lie near the current sheet. This implies that 
high energy particles supporting the pulsar emis- 
sion may be generated not near the regions of high 
Ed but near those with high values of £^|| . 

4. Solutions Vifith J = pc 

As argued above, the time evolution of the 
MHD equations requires a prescription that re- 
lates the current density J and the fields, given 
above by the different forms of Ohm's law. How- 
ever, given that any value of E|| appropriate for 
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pulsars would quickly accelerate any charges to 
velocities close to c, it would be reasonable to as- 
sume that the relation J ~ pc will be attained 
between the current and charge densities. Moti- 
vated by these considerations, rather than deter- 
mine the current J from its relation to Ey, we 
determine it from its relation to the local charge 
density p. Specifically, we focus our attention to 
solutions with J / pc — 1. This specific ratio is 
special because it discriminates between magneto- 
spheres consistent with only a single sign charge 
carrier (J/pc ^ 1) and those that require carri- 
ers of both signs of charges (J/pc > 1) to ensure 
the counter streaming required in space-like cur- 
rent flows. Since each flow line with J/pc > f 
requires at some point a source of charges of the 
opposite sign, it is considered that such flows by 
necessity involve the production of pairs and their 
associated cascades. 

One well-studied location for charge produc- 
tion in the pulsar magnetosphere is the polar cap, 
where particle acceleration and pair cascades are 
thought to take place. If charges of either sign 
are freely supplied by the neutron star surface, 
as is thought t o be the case in all but magnetar- 
strength fields ( Medin fc Lai 2010f) . then a max- 
imum current of J/pc = 1 can flow outward 
above the polar cap in Space-Charge Limited Flow 
(SCLF). In this case, an accelerating electric field 
E\\ is produced by a charge deficit Ap = (p — pgj) 
that develops because the GJ density is differ- 
ent fr om the actual charge flow along the field 
lines ( Arons fc Scharlcmann 197^. This electric 
field then supports a current of primary electrons 
that radiates 7-rays, initiating pair cascades. The 
pairs efficiently screen the above a pair for- 
mation front (PFF), by downward acceleration of 
positrons. This charge density excess is small rela- 
tive to the total charge, i.e. Ap/p = e ^ 1 and the 
current that flows into the magnetosphere above 
the P FF has J/pc ~ 1 ([Harding fc Muslimov 
200 Therefore the steady-state SCLF accel- 



eration models arc only compatible with a small 
range {J/ pc ~ 1) of the cur rent to charge den- 
sity ratio. iBeloborodovl ( 20081 ) demonstrated that 
no discharge occurs if J/pc <1 and that the E\\ 
and the discharge generated when J/pc > 1 or 
J/pc < (implying that the current direction 
is opposite to the one that corresponds to out- 
ward movement of the Goldreich- Julian charges) 



are strongly time dependent. This has been con- 
firmed b y models of time- dependent SCLF pair 
cascades (iTimo khin fc Aro ns 2011). 

With this emphasis on the ratio J/pc = 1, 
we show in Figure [TOl the values of J/pc in color 
scale together with the magnetic field lines on the 
poloidal plane for the parameters of the simula- 
tions of prescription (B) given earlier in Figs. [Si- 
ll] and prescription (C) given in Figs. EHU These 
plots can be directly compared to the FFE mag- 
netosphere of Fig. [1] 

The first two columns of Fig. [10] indicate that 
as cr in prescription (B) ranges from large values 
((7 = 24r2; the FFE has strictly speaking a — 00) 
to (vacuurn) the ratio J/pc decreases tending to 
for cr — > 0|j This can be seen by inspection of 
Eq. ([5]), the first term of which corresponds always 
to a time-like current, with a space-like current 
arising only through the effects of the second term. 
Hence there must be some spatial distribution of 
a values below which no space-like region exists 
because the crE\\ term is small enough. 

The next two columns of Fig. [TOl showing the 
J/pc ratio of Prescription (C) - SEE, are interest- 
ing because this prescription leads to J/pc — 1 
for (T ^> 0, while tending near the FFE values 
for (T — 00, but only in the space-like current re- 
gions of the FFE solutions. The time-like current 
regions of the FFE solutions are now only "effec- 
tively time-like" , as discussed by Gruzinov, result- 
ing in the "fluctuating" or "noisy" behavior that 
is apparent there. 

Suspecting that the "fluctuating" (or "noisy") 
behavior of the effectively time-like current regions 
of SEE is due to the imposed relation between J 
and pc, we modified prescription (B) to search for 
local values of a that would be consistent with 
J/pc ~ 1. More specifically, we consider a a vari- 
able in time and space so that the total J provided 
by Eq. ([S]) is everywhere and at every time step, 
equal to pc. However, we found that there are re- 
gions where goes to zero for a current J less 
than pc. In these regions the second term in the 
expression for J, namely cr , tends to the second 
term of the FFE current expression (Eq. [3]) be- 
ing unable to comply with the J — pc constraint. 
These regions demand extremely high values of 



^The oblique rotators have J/ pc ratios near for mucli lower 
values of cr than those of the first column of Fig. 1101 
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a in order to achieve the J / pc = 1 relation due 
to the vanishing iSy values. We have dealt with 
this numerically by setting a maximum value of a 
{— 50ri) that provided a realistic suppression of 
the in these regions. 

Figure fTTl shows the results of such an approach. 
The first three columns provide respectively the 
values of the poloidal current Jp, charge density p 
and parallel electric component and are simi- 
lar to those presented in Figs. [31 SI [71 El while the 
fourth column shows the values of the J / pc ratio 
in color scale together with the poloidal magnetic 
field lines. The corresponding L and Ed values 
as well as their ratios Ed/L have been denoted 
in Fig. [H by the horizontal dashed lines. Not sur- 
prisingly, these solutions are similar to those of 
prescription (C) - SFE in the cr — ;> regime (see 
Fig. [8]). The essential difference between these 
two approaches is that, in distinction with SFE, 
this approach can handle time-like currents with- 
out the need - like SFE - to do so only on average 
(see blue and purple color regions in the fourth 
column) . 

To allow our solutions even larger flexibility in 
attaining J = pc we have looked for prescriptions 
that divorce J from E||, replacing the second term 
of Eq. ((SJ , with a term of the form 



J|| = ,/ B, 



(12) 



where / is a scalar quantity having a local value 
that ensures the local value of the ratio J/ pc to 
be equal to 1. The difference from the previous 
approach and that of prescription (B) is that this 
term is not related anymore to the value of iSy, 
other than that it is directed along (or against for 
/ < 0) the direction of the magnetic field. The 
sign of / is determined so that the direction of J|| 
coincides with that of E|| . The results of these sim- 
ulations are shown in Fig. [12] the panels of which 
are similar to those of Fig. [TI] The global struc- 
ture in this case is very similar to that of the simu- 
lation presented in Fig. [TT] except for the time- like 
regions (especially those near the star) that now 
exhibit the "noisy" behavior of the SFE simula- 
tions. However, the ratio J/ pc is now everywhere 
equal to 1 . Our solution has therefore the required 
property {J/pc = 1) at the expense of producing 
a "noisy" behavior indicative of time variations at 
these limited spatial scales. 

This oscillatory behavior (both here and in the 



SFE cases) is a rather generic feature that occurs 
in the time- like current regions for values of J/pc 
larger than those needed to "kill" the local E|| im- 
plied by the field evolution equations; both within 
the prescription of Eq. and within SFE, at 
these higher values, J and E^\ end up pointing in 
different directions and an oscillatory behavior en- 
sues. Smaller values of J/pc, cannot drive E|| to 
zero and allow the magnetosphere to attain steady 
state, leaving a non-oscillatory component of E|| 
along the direction of J||. Viewed this way, the 
ideal condition E ■ B = is one that requires a lo- 
cal delicate balance between its parameters (e.g. 
J, p) . If, for whatever reason or necessity, this bal- 
ance breaks, the magnetosphere develops locally 
an oscillatory behavior. The robustness and the 
true evolution of this behavior will be determined 
by the microphysical plasma properties as they 
mesh within the global framework of the magne- 
tosphere electrodynamics. 

The conclusion of this exercise is that solutions 
with a specific value of J, for example J — pc ev- 
erywhere, can be implemented with a sufficiently 
general expression for J|| , provided that one is will- 
ing to accept locally time varying solutions. 

Note that in this prescription, as well as in 
SFE, the current is forced to follow the relation 
J — pc 01 J > pc (for (T ^ in SFE) even in the 
closed field line regions producing even there this 
"noisy" (fiuctuating) behavior. Although apply- 
ing this J — p prescription there is not forbidden a 
priori, it would be reasonable to consider that this 
condition is actually not valid in the closed field 
lines and that a different condition should be used 
there. However, the implementation of a different 
prescription for these regions remains cumbersome 
because it demands that we determine whether the 
magnetic field lines at each grid point are actually 
open or closed and then apply for them this dif- 
ferent condition. Nevertheless, even if we were to 
enforce a prescription that removes the fluctuat- 
ing (noisy) behavior in the closed field line regions, 
this should not affect the similar behavior on the 
open field line regions. 

5. Discussion and Conclusions 

Up to now, the pulsar radiation problem 
has been stud ied considering either the vacuum 
(lDeutschlll955^ or the FEE solutions (jContopoulos et al 
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1999 : Spitkovskvl2006t Kalapotharakos fc Contopoulos 



20091 ) for the structure of the underlying magne- 
tosphere. The vacuum solutions have analytical 
expressions while the FFE solutions have only 
been studied numerically in the last decade. How- 
ever, neither solution is compatible with the ob- 
served radiation; the former because it is devoid 
of particles and the latter because it precludes the 
presence of any E||. The physically acceptable so- 
lutions have to lie somewhere in the middle. As 
of now, there are no known solutions that incor- 
porate self-consistently the global magnetospheric 
structure along with the microphysics of particle 
acceleration and radiation emission. The goal of 
the present work has been to explore the proper- 
ties of pulsar magnetospheres under a variety of 
prescriptions for the macroscopic properties of the 
underlying microphysics in anticipation of a more 
detailed future treatment of these processes. It 
should be noted that these processes lie outside 
the purview of the equations used to evolve the 
electric and magnetic fields. Nonetheless, for most 
models examined the dissipation power does not 
exceed 10-20% of the spin-down one; this value 
is consistent with the observe d radiative efficien - 
cies of the millisecond pulsars ( Abdo et al.ll2009a ). 
that are the ones reaching the highest efficiencies 
of all pulsars, suggesting that energetically, realis- 
tic magnetospheres may be well approximated by 
our models. 

In order to explore the non-Ideal E ■ B 7^ 
regime we considered several prescriptions dis- 
cussed in §3, 4 and analyzed their properties in 
detail (Figures 2-9). Despite their intrinsic in- 
terest, the solutions discussed above still need to 
be examined more closely for consistency with the 
microphysics responsible for the set up and clo- 
sure of the corresponding magnetospheric circuits, 
which are characterized by a considerable range in 
the distribution of the ratio J / pc. On the other 
hand, there are studies arguing that it is essential 
that the ratio J / pc be near 1 in order for steady- 
state pair cascades to be generated near the po- 
lar cap. Besides that, solutions corresponding to 
values J I pc < 1 can be considered as the range of 
configurations without pair production, since such 
solutions can be supported by charge separated 
flows. 

With the value J/pc = 1 given this particular 
significance, we have pursued a search for solu- 



tions that are driven by adherence to this condi- 
tion in as broad a spatial range as possible. As 
noted in the previous section we did so by search- 
ing for solutions of variable conductivity ct, or by 
employing expressions for the parallel current un- 
constrained by the values of the parallel electric 
field. We found that arbitrary values of cr, while 
allowing J j pc — 1 in most areas, in certain others 
can drive E|| to zero sufficiently fast to restrict J 
to values J < pc. As noted in the previous sec- 
tion, a more general expression for J not related 
to E|| can produce J / pc = 1 everywhere but only 
at the expense of fluctuations in regions where the 
value of J/ pc is larger than that necessary to null 
the local E|| implied by the held evolution equa- 
tions. These results indicate the importance of 
the value of the ratio J/ pc vis-a-vis the presence 
of pair cascades and steady-state emission in pul- 
sar magnetospheres. 

The actual behavior of the magnetosphere is 
likely to be more complicated than that described 
by the current density prescriptions provided in 
this work fEas.[51[5]and [T^ . Physics beyond those 
of the magnetic and electric field evolution equa- 
tions will introduce novel time scales that will de- 
termine the response of the current density to the 
parallel electric field component and vice versa. 
The implications of such an oscillatory behavior 
which involves an electric field that alternates di- 
rection, may in fact provide an account of the pul- 
sar coherent radio emission, as such a behavior is 
the essential ingredie nt of cyclotro n maser emis- 
sion ( Melrose 1978; Levinson et al. : 2005; Mclrosi 
20061: ILuo fc Mehosd I2OO8I ). To the best of our 



knowledge, the above argument is the first ever to 
suggest a relation between the global properties 
of pulsar magnetospheres and their coherent radio 
emission. 

These more general models of non-ideal mag- 
netospheres discussed above are possibly related 
to the behavior and properties of intermittent 
pulsars, which may, if more of them are discov- 
ered, provide meaningful constraints on their ra- 
diation mechanism and magnetospheric dynam- 
ics. Intermittent pulsars are pulsars found in "on" 
and "off' radio states ( Lvnd [20091 ) . The inter- 
esting aspect of these objects is that they ex- 
hibit different spin-down rates in their "on" and 
"off' states; even more interesting is the fact 
that they apparently "remember" the spin-down 
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rate they had in the previous cycle, as they re- 
turn to the same spin-down rate when switching 
from "on" to "off" and vice versa. The ratios of 
the "on" to "off" spin-down rates are 1-5 - 1.7 
for pulsars PSR B1931 -H24 (iKramer et a"D '20061) 
and P SR J1832-h0029 (jlorimer et al.l |2006; Ly3 
2009I ). while it is 2. 5 for the more recen tly observed 
PSR J1841-0500 (|Camilo et al.ll2012h . Interest- 



ingly these values are well below the minimum 
ratio 3 between the FFE and the vacu um spin- 
down rates. Based on this observation, iLi et al 



(|2011[ ) argued that the "on" state can be the FFE 
solution while the "off" state simply a resistive 
solution corresponding to a certain value of a. 

An alternative view is to identify the "on" state 
with magnetosperes with J/ pc> 1 and the "off' 
one with those of J/pc < 1, given the necessity 
of the presence of pairs (accompanied by the ra- 
dio emission) in the former and their likely ab- 
sence (and also of radiation) in the latter. In fact, 
the ratio of spin-down rates between the FFE and 
the J/pc ~ 1^ magnetospheres is within 1.4 — 1.8 
(see Figs. [2] and [6|), with the higher values corre- 
sponding to the lower values of a. This would be 
consistent with the observations of the first two in- 
termittent pulsars. However, this is not consistent 
with the spin-down ratio of the third intermittent 
pulsar whose "on" and "off" states may be magne- 
tospheres with different values of a. The discovery 
of a larger sample of intermittent pulsars (espe- 
cially if one could also put some limits on their 
inclination angles) would be extremely useful in 
constraining magnetospheric models. For exam- 
ple if the ratio of the "on/off" spin-down rate was 
found not to exceed the value 3, one would have to 
exclude the vacuum solution being the "off" state, 
since for sufficiently small a this ratio becomes ar- 
bitrarily large. 

The solutions we presented in this work, be- 
sides their intrinsic interest, can serve to produce 
model pul sar 7— ray light curves in a way similar 
to that of Contopoulos fc KalapotharakosI (2010); 
Bai fc Spitkovskvl (l2010al) : IVenter et all |2009); 
Harding et al.l ( 201lh . Such models will provide a 
means of connecting the structure and physics of 
magnetospheres to observations, since they con- 
strain the pulsar magnetospheric geometry as well 
as the regions of anticipated particle acceleration 
and even the values of the accelerating electric 
fields. These will therefore provide independent 



tests of these models for real pulsars. 

The next level of this study should include mi- 
crophysics at a level sufficient to allow some feed- 
back on the global solutions i.e. the production 
of pairs whose effects are included in the com- 
putation of the local E\\ self-consistently. It will 
be of interest to examine whether such solutions 
are consistent with the charge separated flows or 
whether they require an inherently time dependent 
magnetosphere, a feature that could be tested ob- 
servationally. The time dependence should lead to 
oscillations of the pulsar emission with frequency 
of order v > il/27r for emission near the LC and 
V ^ J7/27r for emission near the polar cap. We 
have searched and found not these oscillations in 
the X-ray domain of several pulsar light curves. 
Either they were not present at the objects we 
searched or they lose their coherence at these much 
lower energies. Nonetheless, we intend to study 
the possibility of the presence of such oscillation 
in a future work. 

We would like to thank the anonymous referee 
for constructive comments that helped us to im- 
prove this paper. 
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Fig. 1. — The FFE solutions for the inclination angles a indicated in the figure. Each panel shows the 
poloidal magnetic field lines (white colored lines) and the ratio J j pc in the color scale shown in the figure. 
We note that the color scale saturates below the value 0.5 and above 1.5 so that details around the value 
1 become evident. The dashed vertical lines denote the Light Cylinder (LC). The length unit is equal to 
-Rlc- We observe that the magnetic field lines open beyond the LC. The fraction of the polar cap region 
with space-like currents {J j pc > 1) increases with a so that, while in the aligned rotator almost the entire 
polar cap is filled with time-like currents {J/ pc < 1), in the perpendicular rotator the entire polar cap is 
filled with space- like currents {J/pc > 1). 



13 



0.01 O.Oj 0.1 0.3 1 3 10 30 




oLj I . . -I n 

0.01 0.03 0.1 0.3 1 3 10 30 

alQ 



Fig. 2. — (a) The Poynting fluxes L (measured on 
the surface of the star) as a function of a in log- 
hnear scale for prescription (B). Red, green and 
blue colors correspond to a = 0° (ahgned), a = 45° 
and a = 90° perpendicular rotator, respectively. 
The horizontal solid hne elements denote the L 
values corresponding to the vacuum (lower value) 
and the FFE (higher value) solutions. Note that 
the L value for a = 0° is similar to the vacuum one 
for a < O.Sri (ri is the angular frequency of the 
star) and it reaches that of the FFE solution only 
for much higher a values; (b) The dissipation en- 
ergy rate Ed integrated over the volume bounded 
by radii r\ — r^, = 0.3i?LC and r2 — 2.5i?LC as 
a function of cr; (c) The fraction Ed/ L never ex- 
ceeds the values 10%-20%, while for cr and 
cr — > oo it goes towards 0. The dashed horizontal 
lines in all three panels denote the values corre- 
sponding to simulations with J / pc = 1 (see §4, 
5). 
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Fig. 3. — The non-Ideal magnetospheric solutions in the poloidal plane (/x, Jl) for prescription (B) and for 
a high a value {a = 2Afl). Each row shows the solutions corresponding to the indicated inclination angles 
a. The left-hand column shows the poloidal current modulus (in color scale) together with the streamlines 
of the poloidal current. The middle column shows the charge density (in color scale) together with the 
field lines of the poloidal magnetic field. The color ranges purple-green and green-red indicate negative and 
positive charge density, respectively. The right-hand column shows the parallel electric field component .Eji 
(in color scale) together with the lines of the poloidal magnetic field. The color ranges purple-green and 
green-red indicate antiparallel and parallel directions of (relative to the magnetic field) respectively. Note 
that the color representation for saturates beyond the absolute value 0.2Blc/c where Slc is the value of 
the stellar magnetic dipole field at the distance Rhc- The structure shown in the first two columns is quite 
similar to that of the FFE solutions (the main difference from the FFE solution can be seen in the third 
column) . 
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Fig. 4. — Similar to Fig. |3] but for a much lower a value (cr = l.Sfi). The main magnetospheric features such 
as separatrices and current sheets can stih be observed, even though they now appear much weaker. The 
magnetic field lines close well beyond the LC and the parallel electric field components reach higher values 
than those of Fig. [3] 
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Fig. 5. — Similar to the third columns of Figs. [31 
|3]but for the prescription (A) and for the b = 0.75 
value. For this value of b the corresponding L 
values are the same to those of the simulations 
presented in Fig. |4l The structure of these solu- 
tions seems very similar (qualitatively) to those of 

Fig.m 
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Fig. 6. — Similar to Fig. [2] but for prescription 
(C) - SFE. This prescription does not cover the 
entire spectrum of solutions between the vacuum 
and FFE ones since for cr — we stiU get a dis- 
sipative configuration rather than the vacuum so- 
lution. However, as cr — > oo, the corresponding 
Poynting flux values tend to those of the FFE so- 
lutions. The dissipative energy loss rate Ed ex- 
hibits a maximum value for cr — > and remains 
non-zero even for the high a values due to the per- 
sistent dissipation on the current layer outside the 
LC. This result seems to confirm Gruzinov who 
argued that the SFE solutions for a ^ oo present 
finite dissipation on the current layer. 
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Fig. 7. — Similar to Figs. [21 [3] but for prescription (C) - SFE and for a high a value {<j — 8^1). The Poynting 
flux of these solutions tend to those of the FFE ones. Nevertheless, we still see magnetic field lines closing 
well outside the LC (especially for a = 0°). The electric current in the SFE prescription is everywhere space- 
like {J/ pc > 1). However, there are regions that are effectively time-like (as defined in the text). These are 
the noisy regions in the polar cap vicinity (for a ^ 90°) seen in the third column. Their noisy character is 
due to the continuously changing direction of the electric field. 
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Fig. 8. — Similar to Fig. [7] but for cr = value. The modulus of the parallel electric field components 
increases along with the area of their influence. The closed field line regions of the oblique rotators show 
significant parallel electric field components; however, the dissipation there is insignificant because the value 
of the product J • E is small. 
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Fig. 9. — The third column of Fig. [8] considering 
the average values of the fields within a stellar pe- 
riod. The "noisy" behavior shown in Fig. [8] is now 
disappeared. These regions reveal their effectively 
time-like behavior. 
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Fig. 10. — The ratio J / pc (in color scale) together with the poloidal magnetic field lines for the solutions 
presented in Figs. [3l HI [71 [8] as indicated in the Figure. The ratio J/ pc in prescription (B) decreases in 
general as a decreases and goes towards for a ^ Q. In prescription (C) - SFE the ratio J / pc goes to 1 
(null current) as cr — > while it is higher than 1 for cr > 0. However, in both these cases there are regions 
where the current has on average time-like [J / pc < 1) behavior (effectively time-like; see Fig. [9]). 
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Fig. 11. — The solutions employing prescription (B) with the cr— values adapted locally to achieve 3 j pc as 
close to 1 as possible. The first three columns are similar to those of Figs. |3l IH [71 18] while the fourth column 
shows the corresponding ratio J j pc in color scale. The main features (e.g. current sheets, separatrices) of 
the FFE solutions are still discernible. We observe parallel electric field components across the polar caps 
and along the separatrices. For the oblique rotators parallel electric field components exist also in parts 
of the closed field line regions. However, the latter regions show low dissipation energy rate Eu since the 
corresponding product J • E is small. In the fourth column we observe regions with J / pc < \. In these 
regions a current density J smaller than the value pc is enough so that the corresponding parallel electric 
field component vanishes. This kind of solutions are the upper limit of the solutions that can be considered 
that are supported by charge separated flows. 
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Fig. 12. — Similar to Fig. [11] but with the parahel current density independent of the modulus of E|| (see 
Eq. [T2|) . As shown in the fourth column, this prescription achieves J/pc = 1 in the entire magnetosphere. 
This leads to an oscillating behavior for E|| in the time-like current regions of the simulation presented in 
Fig. [m similar to that of SFE (Fig. |8]). We note that this effect is more general and appears wherever 
the ratio J/pc is enforced to have values above those that are sufRcient to cancel out the parallel electric 
component. This behavior is the result of the mutual response between the electric field and the current 
density and may be crucial for the coherent radio emission. 
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